Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN I^TeX style file v2.2) 



A giant molecular cloud falling through the heart of 
Cygnus A: clues to the triggering of the activity 



o 
o 

(N 
>^ 



> 

(N 
m 

o 

o 

O 



M. J. Bellamy/* and C. N. Tadhunter.^ 

^Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, UK 



To be submitted for publication in the Monthly Notices of the Royal Astronomical Society 

2 February 2008 



ABSTRACT 

We present intermediate resolution near-IR long-slit spectroscopic data for the nearby 
radio galaxy Cygnus A (3C 405) (obtained with the NIRSPEC spectrograph on the 
Keck II telescope). The data reveal considerable complexity in the near-IR emission 
line kinematics, including line splittings of 200-350 km s~^ and a mixture of narrow 
(FWHM -200 km s"^) and broad (FWHM ~700 km s"^) components to the emission 
lines. It is notable that the Paa and H2 emission lines show markedly different kine- 
matics, both on- and off-nucleus. Overall, the data provide evidence for the presence 
of a giant molecular cloud falling through the heart of the Cygnus A host galaxy, the 
motion of which is not driven by the AGN itself. We suggest that this cloud may be 
connected to the triggering of the activity in this highly powerful AGN. We also detect 
split H2 components on the nucleus that are likely to originate in the circum-nuclear 
torus. 

Key words: galaxies: active - galaxies: individual: 3C405 - infrared: galaxies - 
quasars: emission lines - quasars: general 



1 INTRODUCTION 

The prototype radio galaxy Cygnus A {3C 405) is the most 
luminous extragalactic radio source in the local universe 
and the only nearby (z=0.0558) radio galaxy of compara- 
ble power to 3C sources at z~l. In consequence, it is one of 
the most well-studied radio galaxies. Previous work has con- 
centrated on using Cygn us A to test the unified schemes for 
powerful radio gal axies jUeno et alJll994 lOele et al]ll997t 
lYoung et al]l2002l) and to search for signs of the impact of 
the activi ty on the host galaxy in the form of AGN-driven 
outflows iTavlor. Tadhunter fc Robinsonll2003^ . Dynamical 
studies ha ve also led to estimates of the supermassive black 
hole mass (iTadhunter et al.l^2003^ . 

However, despite the large amount of progress that has 
been made over the last decade in the understanding of 
Cygnus A and its host galaxy, relatively little is known 
about the events that triggered the activity. In this paper 
we present near-IR spectroscopic observations, of higher res- 
olution than previously published near-IR da ta (Ward ct al. 
' I99l|: iThornton. Stockton fc Ridewavl Il999l : Fwilman et alJ 
20001) . that provide evidence for non-equilibrium gas mo- 
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tions that may have a direct bearing on the triggering of the 
activity in this poweful AGN. 

In order to put our observations in context, we start by 
summarising our current state of knowledge of the geometry 
and kinematics of the near- nuclear regions of Cygnus A. 



1.1 Geometry and inclination of the near- nuclear 
regions 

HST imaging has revealed an edge-brightened biconi- 
cal structure on t he scal e of a f ew hundred parsecs 
J Jackson et"al]|l99a Tadhunter et aljFl999.1 This structure 
is most likely a direct result of the nuclear-driven winds 
hollowing out 'ionisation cones' in the kpc-scale dust lane 
in which the AGN is embedded. Moreover, the images also 
show evi dence for a starburst ring associated with the dust 
lane tFosburv et ali . 19991. From these data, the best esti- 
mates for the geometry of the biconical structure are open- 
ing half-angle 6i ~ 60° and inclination of the cone axis rel- 

2 

ative to the line-of-sight i ~ 30°with the NW cone oriented 
towards us and the S E cone oriented away from us (see 
ITadhunter et al]l2003l for a full discussion). From this ge- 
ometry we know that we are likely to be viewing Cygnus A 
very close to the opening angle of the NW cone, but we also 
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know that i < (90 - di) as we do not have a direct view of 
the quasar nucleus at optical wavelengths. 

1.2 Kinematic components in the near-nuclear 
regions 

The central regions of Cygnus A show considerable com- 
plexity in their emission line kinematics. The various kine- 
matic components can be tied to different processes associ- 
ated with the AGN as follows: 

Gravitational motions 

Kpc-scale disc: near-IR Keck II spectroscopic data shows 
evidence for a rotating disc associated with the 3 kpc scale 
dust lane. The gas is in Keplerian rotation about a stellar 
core and an unresolved point mass of (2.5 ± 0.7) X 10^ solar 
masses. The rotation axis of the gas is aligned at ~9 °to that 
of the large-scale radio axis jTadhunter et alJl2003l) . 

300pc-radius disc: Higher spa t ial re solution HST/STIS 
optical data of iTadhunter et alJ i2003l) show evidence for 
a smaller disc structure. The lines from this disc are sig- 
nificantly broadened throughout (FWHM 500-900 km s"^). 
Although the slit positions closest to the nucleus show ev- 
idence for rotation about a central black hole with mass 
similar to that deduced from the Keckll/NIRSPEC data, 
there is evidence for deviations from pure circular motions 
in the NW cone. 

Evidence for outflovirs 

Emission line outflows: Intermediate spe ctral resolu- 
tion optical data of the [O IIIJA5007 line (iTavlor et alJ 
|2003^ provide evidence for outflowing gas in the ionisation 
cones. A ~300 km s~^ outflow is detected in the NW cone 
(blueshifted), whereas a ~400 km s~^ outflow is detected in 
the SE cone (redshifted) . These results are consisitent with 
gas being driven out of the cones by AGN-induced winds. 

Extreme outflow: an extreme outflow component is de- 
tected in [O III1AA(5007,4949) emission in the NW cone 
jTadhunteil 179911. This component lies on the radio axis 
itself and is blueshifted by 1300 - 1800 km s~^; it is likely 
to represent gas entrained in the outflowing radio plasma. 

Scattering outflow: The redshifted [O III]A5007 feature 
de tected in polarised light, in both cones and the nucleus, 
bv lvan Bemmel et alJ i2003l) is consistent with the presence 
of a scattering outflow with velocity in the range 150-450 
km s~^. 

Evidence for inflows 

Infalling H I: Two components of H I 21cm absorp- 
tion are detected in t he nuclear regions of Cygnus A 
JConwav fc Blancdll995ll . In the rest frame, defined by the 
narrow component of Paa in the nuclear aperture (see be- 
low), these components have redshifts of 227 ± 9 km s~^ 
and 48 ± 9 km s~^. Since these absorption components are 
detected against the radio core, and are therefore in the 
foreground, they must be associated with infalling material. 



Although there is now substantial evidence for the gravi- 
tational motions and AGN-induced outflows in Cygnus A, 
the evidence for inflowing material is re latively sparse and 
rests solely in the H I observations of IConwav fc Blanca 
(1995). However, it is the inflow component that is most 
likely to provide the strongest clues to the triggering of the 
activity. Therefore it is important to obtain further informa- 
tion about the kinematics and distribution of t he inflowing 
gas. In this context we note that lCanalizo et al.l (2003) have 
recently reported the discovery of a red, secondary point 
source in the nuclear regions, ~400pc southwest of the radio 
nucleus. They argue that this secondary nucleus represents 
the debris of a minor merger that may be fuelling the AGN 
activity. 

We assume the cosmological parameters of Ho =75 km 
s~^ Mpc~^ and qo=0 throughout this paper. For these pa- 
rameters 1.00 arcsec corresponds to 1.00 kpc at the redshift 
of Cygnus A. 



2 DATA COLLECTION AND REDUCTION 

The K-band near-IR spectra of Cygnus A were taken on 2000 
May 22 using the NIRSPEC spectrograph in grating mode 
on the Keck II telescope at the Mauna Kea complex. Long- 
slit spectra were obtained with the slit oriented along the ra- 
dio axis (PA105) and along PA180: in each case the slit was 
centred on the nucleus. Velocity curves derived from single- 
Gaussian fits to these data - which provi de evidence for a 
rotating kpc-scale disc - are presented in ITadhunter et al] 
(2003). In this paper we will be concerned only with the ra- 
dio axis position angle (PA105). Four exposures of 400 sec- 
onds were obtained (at each slit position) with the exposures 
nodded in an ABBA pattern to facilitate sky subtraction. 

The spectra were rectified using purpose-written iraf 
routines, and wavelength-calibrated in iraf using an ar- 
gon/neon arc lamp exposure taken at the time of the ob- 
servations. Following the manual removal of cosmic rays us- 
ing the CLEAN routine in FIGARO, the separate exposures 
for each PA were combined to produce a sky-subtracted 
two-dimensional spectrum. Telluric features were removed 
by dividing a high S/N spectrum of the AO standard star 
HD203856, which was observed at the same airmass, into 
the galaxy spectra. Use of the 0.57-arcsec slit resulted in a 
spectral resolution of 10.9 ± 0.5 A (140-170 km s"^ for these 
data), with the observations covering the wavelength range 
19100-23000A. From measurements of stars along the slit, 
we estimate a plate scale of 0.178 ± 0.004 arcsec prxel"^, 
and an effective seeing of 0.75 ± 0.05 arcsec (FWHM) for 
the observations. Measurements of night-sky emission lines 
in the reduced frames demonstrate that the wavelength scale 
is accurate to within ± 0.5 A (± 8 km s~^) along the full 
length of the slit. 

The exposures of HD203856 were used to fiux-calibrate 
the data, with the assumption that the intrinsic spectral 
energy distribution (SED) of the star is that of a perfect 
black-body at T = 9480 K. The m agnitude-to-flux conver- 
sion was performed with reference to lBessell. Castelli fc Ple'3 
(1998). The Star link FIGARO and dipso packages were used 
to analyse/reduce the data and fit the spectral lines. 
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Table 1. Results of the DIPSO fits to tiie spectral lines in the PA105 central aperture. The FWHMs have been corrected for the 
instrumental width. The redshifts have been heliocentrically corrected. 
"Fit from the H2 !^=l-0 S(l) components. 

'Fit from the Fact components^ 

'=Rest frame wavelengths from lThornton et al.l ll999t) . 
''Rest wavelength not known accurately from atomic physics. 

IWilman et aII(l2000D . the Paa and molecular hydrogen lines 
exhibit different properties and are therefore discussed in 
separate sections. Where velocity scales and shifts are given, 
the zero velocity corresponds to the heliocentric z=0. 05583, 
the centre of the Paa narrow component in the nuclear 
aperture. This systemic redshift is consistent with, though 
slightly lower than that deduced from the PaQ:/H2 rotation 
curves (iTadhunter et al.l 2003: z=0. 05592 ± 0.00005). 



3.1 Paschen alpha 

The Paa emission exhibits kinematics, both on- and off- 
nucle us, that show simila rities to those of the optical [O III] 
lines jTavlor et alJl2003^ . Unlike the case of the molecular 
hydrogen lines (see section 3.2), there are no comparable 
recombination lines with which to check the accuracy of the 
component fittings. However, the Paa line is the brightest 
line in the spectrum, and has a high S/N ratio, so we can 
be confident in the results of the fittings. 

In the central nuclear aperture, the Paa line clearly 
has broad wings that cannot be fitted with a single Gaus- 
sian profile. Fitting shows a central narrow component with 
FWHM — 250 km s~^ and a much broader component with 
FWHM = 750 km s'^ . HST/STIS data of the nuclear re- 
gions shows a rotation-curve in the [O HI] line. Integrated 



3 RESULTS 

To examine the kinematics along the slit with good S/N, 
a number of apertures were extracted. Two of these aper- 
tures reveal the complex kinematics that are discussed in 
this paper. The first aperture is centred on the nucleus (to 
the nearest pixel) and has a width of three pixels (~540 pc). 
The second aperture is centred ~1 880 pc to the NW of the 
nucleus and is four pixels (~720 pc) in width. Fig. 0shows 
the greyscale 2D profiles of the Paa and the H2 z^=l-0 S(l) 
line emission with these apertures overlaid. Fig.|21shows the 
ID spectrum extracted from the PA105 central aperture. 

The emission lines in the aperture spectra were fitted 
with Gaussian profiles using the dipso package. Each line 
was fitted with the minimum number of Gaussians required 
to give an acceptable fit. Table shows the results of the 
line fitting. Several lines show evidence for multiple kine- 
matic components (notably the brighter H2 lines) and these 
components have been fitted separat ely. We see no evidenc e 
for the H2 u=2-l S(2) line seen bv IWilman et al.l ll200Cl) . 
but note that their slit is approximately twice the width of 
the slit used in the acquisition of these Keck II data. This 
may suggest that the H2 1^=2-1 S(2) emission originates in 
a more extended region. 

The Paa and molecular hydrogen lines are discussed 
in more detail in the following. As previously noted by 
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Figure 1. The 2D profiles of the Pa« (top) and H2 u=l-0 S(l) 
(bottom) emission in the PA105 data. The spatial scale runs ver- 
tically with the NW side at the bottom; the velocity scale runs 
horizontally with positive velocities to the right. The extraction 
apertures are marked on the right hand side: a) the nuclear aper- 
ture and b) the extended NW aperture. 



over the inner 600 pc, and fitted with a single-Gaussian, 
the pro file of the [O III1 line gives a FWHM of 850 ± 50 
km s" "^ iTadhunter et al.ll20ol) . The broad Paa component 
is therefore likely to be due to a combination of unresolved 
rotation and turbulent motion in the <300 pc-scale nuclear 
disc. The centre of the broad component is coincident with 
that of the narrow component within la uncertainty. 

The situation is different in the extended NW aperture 
in the PA105 data: this aperture is centred 1 880 pc NW of 
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Figure 2. K-band ID spectrum extracted from the central aper- 
ture of the PA105 data. 
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Figure 3. (Top) The two components of Paa in the extended 
NW aperture (1 880 pc from the nucleus). Marked on the graph 
are the velocity shif ts of the [O III] components in the same region 
iTavlor et al.l^2003^ . (Bottom) The two Pa a components shown 
fitted to the H2 u=l-0 S(l) line in the NW aperture. The relative 
fluxes of the Paa fit are retained as the fitting could not be made 
with these parameters free. 
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the nucleus, along the radio axis in the NW ionisation cone. 
In this aperture the Paa is split into two narrow compo- 
nents (FWHM ~360 km s~^), one of which is significantly 
blueshifted by 240 ± 30 km s~^ in the galaxy-frame (see Fig. 
I^l; the other, brighter, component is redshifted by 120 ± 13 
km s~^. As can be seen in Fig. E] these components provide 
a poor fit to the H2 emission in the same aperture. In the 
same spatial region, two components are seen in the [OIII] 
A5007 line at optical wavelengths: these are shi fted by (-270 
±20) and (-1-100 ± 15) km s"^ respectively llTavlor et alJ 
|20o3)- Allowing for the different slit widths and systematic 
uncertainties, the Paa compo nents are entirely c onsistent 
with the [O III] components of iTavlor et"ai] J2003l) . There- 
fore the blueshifted component is likely to be associated with 
the emission line outfiow in the NW cone. 

3.2 Molecular hydrogen 

The molecular hydrogen emission lines show a complex 
structure of kinematic components that is not seen in any of 
the other emission lines. In the central aperture of the PA105 
data, three distinct components are seen in the H2 lines — 
a central 'narrow' component (FWHM = 160 km s~^) and 
two broader components (FWHMs — 270 and 370 km s~^). 
Fig. 2] shows these 3 components as fitted to the H2 z^=l- 
S(l) line. One of the broader components is blueshifted 
by (370 ± 40) km s~^ with respect to the systemic velocity 
and the other is redshifted by (110 ± 40) km s^^. These 
components were derived from the fitting to the H2 1^=1-0 
S(l) line but also describe very well the profile of the H2 
i'=l-0 S(2) (see Fig. |3J| and, with different intensity ratios, 
the H2 1^=1-0 S(4). They also appear to fit the H2 1^=1-0 
S(3), but this is less convincing due to the blend with the 
[Si VI] line. However, if the H2 v=l-0 S(3) line is fit with 
the H2 i'=l-0 S(l) components, and the [Si VI] line with 
the Paa components, then the blend can be fitted accept- 
ably. Fig. |1] shows that these H2 components do not pro- 
vide a good fit to the Paa emissi on in the central ape rture. 
These results contradict those of lWilman et al.l i2000fl who 
found the ortho-hydrogen lines [S(l) and S(3)] to be sig- 
nificantly broader than the para-hydrogen lines [S(2) and 
S(4)] on the nucleus. This dis crepancy may be a re sult of 
the single- Gaussian fittings of IWilman et all i2000ll being 
insufficient to fit the actual line profiles. 

As can be seen in Fig. Q the molecular hydrogen shows 
a bright, shifted component to the NW of the nucleus. The 
slit is closely aligned with the rotation axis of the kpc-scale 
disc for these observations and, for this reason, the rotation 
of this disc would not lead to such shifted components. The 
extended structure can be seen in all the molecular hydro- 
gen lines in the original 2D spectrum, but the best fit can be 
made with the bright H2 v—l-O S(l) line. Fitting to this line 
on a row-by-row basis reveals two components, one at the 
systemic velocity and one that is significantly redshifted. 
The two components exhibit no significant velocity varia- 
tions across the range of the fitting. Fig.j^shows the velocity 
structure of these components over the 0.8-2 kpc region. Fit- 
ting on a row-by-row basis at projected distances <0.8 kpc 
produces indefinite results without added constraints and 
these fits are therefore not plotted. These results indicate 
the existence of a molecular hydrogen cloud at a projected 
distance of ~1.35 kpc from the nucleus (where the fiux of 
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Figure 4. The three components (see text) of the H2 lines as 
shown fitted to (top) the H2 i'=l-0 S(l) line, (middle) the H2 
1^=1-0 S(2) line and (bottom) the Pao line (all in the nuclear 
aperture). The weak feature redward of the H2 !^=l-0 S(2) line 
is the H2 u=8-6 0(3) line; the weak feature blueward of the Paa 
line is the He I (4F-3D, 3Fo-3D) line. It is apparent that the H2 
components do not provide a good fit to the Paa emission, even 
with the relative fluxes of the components allowed to vary freely. 
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Figure 5. Results of the row-by-row fitting to the H2 v=l-0 S(l) 
hne. Two components are seen in the hne in the NW aperture 
and the variations in the velocities of these components are shown 
between 0.8 and 2 kpc. The points at r = show the velocities of 
the three components detected in the central nuclear aperture. 

the redshifted component peaks), redshifted by 243 ± 7 km 
in the rest frame. This is consistent with the redshift 

of 227 ± 9 km s"^ exhibited by one of the on-nucleus H 

I absorption components seen bv lConwav fc Blancol il995l) . 

The emission also appears to be spatially resolved with the 

flux distribution having a FWHM of ~500pc (assuming a 

Gaussian intensity profile). 

Overall, the H2 emssion shows kinematic components 

that are not seen in other lines. Possible explanations are 

discussed below. 



4 DISCUSSION 

4.1 The nuclear emission 

The nuclear kinematics of the Paa line are undramatic and 
show no unusual or unexpected features. Broad (FWHM = 
750 km s~^) and narrow (FWHM = 250 km s~^) compo- 
nents are detected, each symmetrically centred at the sys- 
temic redshift, with no evidence for any other components. 
T hese line widths are typical of those seen in the optical data 
of lTavlor et al.l ll2003h . The broad component compares well 
with the [O HI] width integrated over the inner 30 pc (ra- 
dius) disc as seen in the HST/STIS data of Tadhunt eretai] 
i2003r , suggesting it represents unresolved rotation and tur- 
bulent motions close to the nucleus. 

The narrow component of the H2 lines, although nar- 
rower than the narrow Paa component (FWHM = 160 km 
s~^), is centred at the same systemic redshift within the un- 
certainties. The red and blueshifted components, however, 
are not so easily understood. Their widths, although broader 
than the central H2 component, are certainly not individu- 
ally comparable to the broad Paa component. Furthermore, 
their velocity shifts are significantly asymmetric about the 
centre. What is the nature of these components? Are we 
seeing outflows in the nuclear regions? This is certainly a 
possibility, but it is hard to explain how we would then 
see the redshifted component. The blueshifted component 



would obviously correspond to material flowing out into the 
cone that is angled towards us and we would have a rela- 
tively clear view of this. The redshifted component, however, 
would correspond to an outflow on the far side of the nucleus 
and, as it is seen in the nuclear aperture, it would necessar- 
ily lay close to the circum- nuclear torus. It would therefore 
be subject to a high degree of extinction, even at near-IR 
wavelengths, and should not be so apparent. The fact that 
the redshifted component is the most luminous of the three, 
and a factor of 5x brighter than the blueshifted component, 
only further compounds the problem. Another considera- 
tion is that, if these components represent outflows within 
the ionisation cones, similar components should be seen in 
Paa due to the ionisation of the outflowing material. Such 
components are not seen. 

It is more likely that these velocity-shifted components 
originate in the circumnuclear torus and are linked to its 
rotation. If, as is believed, we are viewing Cygnus A at close 
to the opening angle of the NW ionisation cone, then there 
is likely to be a relatively clear view of the outer parts of 
the circum-nuclear torus. The question to ask is why do we 
see three distinctly separate kinematic components rather 
than a continuous broadening of the molecular lines due to 
unresolved rotation? It is possible that clumpy obscuring 
material partially blocks our view and so we only sample 
part of the velocity structure. There is good evidence that 
there is patchy dust obscuration within the central region 
JShaw fc T adhuntcr"l994'. Stockt on. Ridewav fc Liihlll994l. 
I Jackson et al. 1996') - although on a larger scale - so this is 
by no means an unreasonable explanation. 

The H2 u=l-0 S(3)/H2 S(l) flux ratio is use- 

ful as a diagnostic indicator of the excitation mechanism of 
the H2 gas. In the nuclear aperture this ratio is found to 
be 1.23 ± 0.03 (found by integrating the flux of the two 
lines over equivalent ra nges). X-ray heating is expected to 
lead to a ratio of 1.2 jLepp fc McCravlll983ll . Our result 
is in good agreement with this and is consistent with the 
nuclear H2 gas being excited by emission from the AGN. 
iThornton et alJ il999l) previously obtained a value of 1.20 
± 0.14 for this ratio; although consistent with our result, 
this value fails to rule out shock excitation, expecte d to re- 
sult in a ratio of 0.8-1.1 jKwani ri977: HoUenbach fc ShulJ 
^73). The more precise value we have obtained seems to 
rule out shock excitation as a plausible mechanism. 

As noted above, the split components visible in the H2 
emission are conspicuously absent in Paa (see Fig. |1J . This 
is explainable, however, by noting that the Paa emission 
from the torus would originate close to the ionised inner- 
edge in a relatively confined region, whereas the molecular 
hydrogen would originate in a larger, more extended region 
beyond this. Because of the large viewing angle with respect 
to the torus axis (^^60°), the nearside inner-edge of the torus 
lies behind the optically thick outer regions and we would 
therefore expect the Paa to be more highly obscured than 
the H2. The Paa from the far inner-edge should be visible 
however, but this will sample only a small velocity range 
and may therefore be detected as part of the observed broad 
component. 
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4.2 The extended H2 emission 

What is the origin of the redshifted, extra- nuclear H2 com- 
ponent as seen to the NW of the nucleus (see Fig. and 
Fig. IHJ? There are four possible explanations that we will 
discuss: a gas outflow from the nucleus; nuclear light scat- 
tered off material further out; stellar excitation of H2 in a 
star-forming region; and an infalling or kinematically decou- 
pled gas cloud passing through or close to the nucleus. 

First we shall consider the possibility that this is emis- 
sion from an outflowing gas cloud. On the face of it, this 
is the most apparent and satisfying explanation: the nu- 
clear winds and radio jets provide ready-made driving mech- 
anisms for outflows. Such outflows are an expected by- 
product of AGN activity and evidence for this in Cygnus 
A has already been detected in the [O III] lines (e.g. 
iTavlor et all l200.?l : Iron Bemmel et alJ l200.t l. Indirect evi- 
dence for this phenomenon in Cygnus A is also provided 
by the hollo wed out cones apparent in near-IR HST images 
of the core jTadhunter et aljll999D . However, there are a 
number of problems with this explanation for the shifted 
H2 component. Firstly, if this cloud is outflowing within the 
ionisation cone, Paa emission from the photoionised edge of 
the cloud facing the nucleus would be observed with similar 
velocity shifts, but this is not detected. If this component 
can be seen in the H2 line it should certainly be seen in 
the much brighter Paa line. Secondly, since the NW cone 
is angled towards us, the far side of the cone lies close to 
the plane of the sky (<30°): any outflowing cloud showing 
a redshift would therefore have to lie close to the far side 
of the cone and the necessary de-projected outflow velocity 
would be in excess of 450 km s~^. It is hard to envisage how 
the gas could be accelerated to such high velocities without 
being shocked, and shock speeds in excess of ~50 km s^^ 
woul d lead to dissociation of the H2 a s well as to Paa emis- 
sion (lKwanlll977l : iDraine et al.l ll98^. The combination of 
these factors renders an outflow as an unlikely explanation. 

There is also the possibility that the extended emission 
is nuclear light scattered off outflowing dust in the cone. In 
this case the geometry will always result in a net redshift of 
the emission unless the outflowin g dust is moving directl y 
towards us along our line-of-sight (Ivan Bemmel et alll2003ll . 
This would explain the redshift, but still does not explain the 
lack of an equivalent Pact component since, to reflect near- 
IR light, the scattering material must be optically thick and 
therefore relatively grey in its spectral response. This expla- 
nation was proposed in ITadhunter et al.1 ll2003l . see section 
5.7) to explain the redshifted [O HI] emission seen in the 
NW cone in HST-STI S optical data; t he spectropolarime- 
try results of^yan Be mmel et al.l i2003r) are consistent with 
this model. However, this [O III] emission lies much closer 
to the nucleus than the H2 emission (~200 pc as opposed 
to ~1 350 pc) and shows a significantly larger velocity shift 
(^^450 km s~^ a s oppo sed to ~200 km s~^). As noted in 
ITadhunter et all ll2003i) . there is also a problem with the 
stability of any such dust within the ionisation cone and 
whether it could exist there without being destroyed, par- 
ticularly since any gas associated with the dust must be hot 
enough to be in an ionised state to avoid being detected in 
emission itself. 

Alternatively, the H2 emission may arise from the pho- 
todissociation zones associated with H II regions in which 



stars are forming; the stars themselves may remain unseen 
inside the optically thick clouds. However, near-IR observa- 
tions of H II galaxies demonstrate that the Paa emission 
from such regions is strong relati ve to the H2, with a ra- 
tio of PaQ/H2 1^=1-0 S(l) ~5-30 JVeilleux et alJll99g^ . In 
contrast, the maximum ratio allowed by attempted three 
component fits to Paa is Paa/H2 i/=l-0 S(l) ~2 for the 
redshifted component of H2. Therefore we conclude that the 
redshifted H2 component is unlikely to be excited by a young 
star forming region. 

Finally we consider the case of the H2 emission be- 
ing associated with an infalling cloud. Indeed, the possi- 
bility that Cygnus A is undergoing a minor merger event 
JCanalizo et all2003f) suggests that such components should 
be present and natuarally explains the fact that the compo- 
nent we see is significantly redshifted. Infalling material has 
al ready been seen in the H I 21 cm absorption components 
of lConwav fc Blancd (|^9^ and the redshift of one of these 
components is in good agreement with the redshift of the ex- 
tended H2 component, suggesting a common origin for the 
two. However, the problem of the missing Paa component 
remains. If it is falling through the ionisation cone, the cloud 
should be photoionised on its AGN-facing side regardless of 
the cloud's origin. Alternatively, if it is outside the ionisa- 
tion cone, and therefore shielded from the AGN, then what 
is stimulating the molecular emission? An intriguing possi- 
bility is that the cloud may be in the foreground outside the 
ionisation cone. Here it would be shielded from the major- 
ity of the AGN emission but would still be illuminated by 
the more penetrating hard X-ray radiation. The cloud will 
appear optically thin to these hard X-rays resulting in less 
photoionisation (and hence little emission from recombina- 
tion lines such as Paa), but with the fast electrons released 
in the rare photoionisation events capable of thermally ex- 
citing large quantities of the molecular gas (LeuD & McCravl 
1983). This is a plausible explanation for the discrepancy be- 
tween the Paa and H2 emission profiles and one that rules 
out a simple outfiow, since it is hard to see how the AGN 
could be driving outfiows beyond the confines of its ionisa- 
tion cones. 

As a test of the infalling cloud premise, the H2 lumi- 
nosity of the cloud can be compared to the hard X-ray lu- 
minosity of the AGN. Approximately 0.2% of the 2-10 keV 
X-ray flux incident on t he gas cloud will be re -emitted in 
the H2 u=l-0 S(l) line jLepp fc McCravl Il983^ . therefore 
the 2-10keV/H2 1^=1-0 S(l) luminosity ratio leads directly 
to an estimate of the covering factor of the cloud. For our 
cosmology, the intrinsic 2-10 keV hard X-ra y luminosity of 
the AGN (deduced from lYoung et al.ll2002^ is 1.64 x lO'*" 
erg s~^. If the molecular cloud is outside the ionisation cone 
it should see this X-ray emission through a similar amount 
of obscuration as we do along our line-of-sight. We therefore 
reapply Young et al.'s inferred extinction (column density 
Nh ~ 2.0 X 10^^ cm~^) to derive an obscured luminosity 
of 1.44 X 10*'' erg s~'. The integrated luminosity of the ex- 
tended H2 1^=1-0 S(l) emission is (4.9 ± 0.3) x lO^'' erg s"\ 
Accounting for the 0.2% efficiency factor, this implies a cloud 
covering factor of ~1.7 x 10~^. At a projected distance of 
1.35 kpc this equates to a cloud with a circular cross-section 
~360 pc in radius. This is large for a typical single CMC, but 
a merger event would be expected to involve a large mass of 
material distributed in several molecular cloud complexes. 
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The extended H2 structure appears to be spatially resolved 
and has an approximate FWHM of -^SOOpc wich is in good 
agreement with the results of the above calculation. 

Certainly, the infalling cloud model seems to be the 
most appealing explanation of the four possibilities. The ob- 
served redshift is readily explained and the lack of a match- 
ing Paa component is not problematic. The missing Paa 
alone is good grounds for eliminating the other explanations 
and, as described, there are also other problems associated 
with them. The fact that the H2 redshift is in agreement with 
the redshift of the H I absorption component is also encour- 
aging. This infalling cloud scenario also tie s in we ll with the 
postulated merger event of ICanalizo et all (|200^. Unfortu- 
nately the S/N of the data is insufficient for the extended 
component to be fit to anything other than the H2 v—l- 
S(l) line, meaning that no information can be obtained 
about the excitation of this extended component. However, 
the motions of the H2 cloud clearly cannot be accomodated 
within the framework of the (equilibrium) kinematics of the 
kpc-scale disc associated with the dust lane. 

Although the infalling H2 cloud may not be responsible 
for fuelling the nuclear activity directly, the process of cap- 
ture of the cloud, and subsequent dissipative settling into the 
kpc-scale disc, will lead to radial gas motions in the disc and 
an increased fuelling rate of the AGN. In such a scenario, 
the dynamical timescale required for all the gas associated 
with the merger event to settle into the disc is likely to be 
comparable with the lifetime of the quasar and jet activity. 



5 CONCLUSIONS AND FURTHER WORK 

Our near-IR data have further strengthened the case for 
the existence of complex, non-gravitational kinematics in 
the core of the nearby radio galaxy Cygnus A. We show 
several distinct Paa components, both on- and off-nucleus, 
that are consi stent with the [O II I1 components seen in the 
optical data of lTavlor et al.l l|20o3)- The molecular hydrogen 
emission, however, shows markedly different kinematics that 
are inconsistent with both the Paa and optical data. 

We have found good evidence for the existence of an 
infalling molecular cloud to the NW of the nucleus and have 
argued that this lies outside the main ionisation cone and is 
being excited by the more penetrating hard X-ray emission 
from the AGN. The redshift of this infalling component is 
in excelle nt agreement with that of the H I absorption de- 
tected bv lConwav fc Blancol (^9^ , suggesting the two phe- 
nomena are linked. The fact that the H I absorption is seen 
on-nucleus, whereas the H2 emission is seen at a projected 
radius of ~1.35 kpc, suggests that the infalling material is 
likely to have a complex and extended spatial structure. 

We suggest the possibility that the non-circular motions 
and infalling material in the nuclear regions are linked to a 
possible merger event iCanalizo et aljf2003f) that may have 
triggered the AGN activity itself. 

However, the problem remains that these long-slit 
data give kinematic information only along the radio axis, 
whereas a minor merger event could be expected to lead to 
non-gravitational motion and flows throughout the cones. 
Future near-IR integral field data that map the structure 
and kinematics of the H2 emission throughout the cone will 
reveal if the redhifted emission is associated with the radio 



jet itself or is more evenly distributed in the near-nuclear re- 
gions. Near-IR adaptive optics spectroscopy of the nuclear 
regions will also be able to clarify the nature of the asymmet- 
rically shifted H2 components seen in our nuclear aperture. 



ACKNOWLEDGEMENTS 

We acknowledge the facilities and thank the staff at the Keck 
II telescope and the Mauna Kea facility. We also thank the 
referee, Neal Jackson, for his helpful comments and sug- 
gestions. MJB is supported by a PPARC studentship. We 
acknowledge the data analysis facilites at Sheffield provided 
by the Starlink Project, which is run by CCLRC on behalf 
of PPARC. 



REFERENCES 

Bessell M. S., Castelh F., Plez B., 1998, A&A, 333, 231 
Canalizo G., Max C, Whysong D., Antonucci R., Dahm 

S. E., 2003, ApJ, 597, 823 
Conway J. E., Blanco P. R., 1995, ApJ, 449, L131 
Draine B. T., Roberge W. G., Dalgarno A., 1983, ApJ, 264, 

485 

Fosbury R. A. E., Vernet J., Villar-Martin M., Cohen 
M. H., Ogle P. M., Tran H. D., 1999, in Rottgering H. 
J. A., Best P. N., Lehnert M. D., eds. The Most Distant 
Radio Galaxies, Proceedings of the colloquium, Amster- 
dam, 15-17 October 1997, Royal Netherlands Academy of 
Arts and Sciences p. 311 
Hollenbach D. J., ShuU J. M., 1977, ApJ, 216, 419 
Jackson N., Tadhunter C, Sparks W. B., 1998, MNRAS, 
301, 131 

Jackson N., Tadhunter C, Sparks W. B., Miley G. K., Mac- 

chetto F., 1996, A&A, 307, L29 
Kwan J., 1977, ApJ, 216, 713 
Lepp S., McCray R., 1983, ApJ, 269, 560 
Ogle P. M., Cohen M. H., Miller J. S., Tran H. D., Fosbury 

R. A. E., Goodrich R. W., 1997, ApJ, 482, L37 
Shaw M., Tadhunter C, 1994, MNRAS, 267, 589 
Stockton A., Ridgway S. E., Lilly S., 1994, AJ, 108, 414 
Tadhunter C, Marconi A., Axon D., Wills K., Robinson 

T. G., Jackson N., 2003, MNRAS, 342, 861 
Tadhunter C. N., 1991, MNRAS, 251, 46P 
Tadhunter C. N., Packham C, Axon D. J., Jackson N. J., 

Hough J. H., Robinson A., Young S., Sparks W., 1999, 

ApJ, 512, L91 

Taylor M. D., Tadhunter C. N., Robinson T. G., 2003, MN- 
RAS, 342, 995 

Thornton R. J., Stockton A., Ridgway S., 1999, AJ, 118, 
1461 

Ueno S., Koyama K., Nishida M., Yamauchi S., Ward M. J., 

1994, ApJ, 431, LI 
van Bemmel I. M., Vernet J., Fosbury R. A. E., Lamers H. 

J. G. L. M., 2003, MNRAS, 345, L13 
Veilleux S., Sanders D. B., Kim D.-C, 1999, ApJ, 522, 139 
Ward M. J., Blanco P. R., Wilson A. S., Nishida M., 1991, 

AJ, 382, 115 

Wilman R. J., Edge A. E., Johnstone R. M., Crawford 
C. S., Fabian A. C, 2000, MNRAS, 318, 1232 



A giant molecular cloud falling through the heart of Cygnus A 9 



Young A. J., Wilson A. S., Terashima Y., Arnaud K. A., 
Smith D. A., 2002, ApJ, 564, 176 



